Synchronization of fractional order chaotic systems * 
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The chaotic dynamics of fractional order systems begin to attract much attentions in recent years. 
In this brief report, we study the master-slave synchronization of fractional order chaotic systems. 
It is shown that fractional order chaotic systems can also be synchronized. 
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Fractional calculus is a 300-year-old topic. Although 
it has a long mathematical history, the applications of 
fractional calculus to physics and engineering are just a 
recent focus of interest [1, 2], Many systems are known 
to display fractional order dynamics, such as viscoelas- 
tic systems [3-5], dielectric polarization [6], electrode- 
electrolyte polarization [7] and electromagnetic waves [8] . 
More recently, many authors begin to investigate the 
chaotic dynamics of fractional order dynamical systems 
[9-17]. In [9], it has been shown that the fractional or- 
der Chua's system of order as low as 2.7 can produce a 
chaotic attractor. In [10], it has been shown that nonau- 
tonomous Duffing systems of order less than 2 can still 
behave in a chaotic manner. In [11], chaotic behaviors of 
the fractional order "jerk" model was studied, in which 
chaotic attractor was obtained with system orders as low 
as 2.1, and in [12] the chaos control of this fractional order 
chaotic system was reported. In [13], chaotic behavior of 
the fractional order Lorenz system was studied, but un- 
fortunately, the results presented in this paper are not 
correct. In [14] and [15], bifurcation and chaotic dynam- 
ics of the fractional order cellular neural networks were 
studied. In [16], chaos and hyperchaos in the fractional 
order Rossler equations were studied, in which we showed 
that chaos can exists in the fractional order Rossler equa- 
tion with order as low as 2.4, and hyperchaos exists in the 
fractional order Rossler hyperchaos equation with order 
as low as 3.8. In [17], we have studied the chaotic behav- 
ior and its control in the fractional order Chen system. In 
[18] , the author present a broad review of existing models 
of fractional kinetics and their connection to dynamical 
models, phase space topology, and other characteristics 
of chaos. 

On the other hand, synchronization of chaotic systems 
has attracted much attentions [19] since the seminal pa- 
per by Pecora and Carroll [20]. In this brief report, we 
study the synchronization of fractional order chaotic sys- 



tems. The analysis of fractional order systems is by no 
means trivial. So, we will numerically investigate this 
topic here. 

There are many definitions of fractional derivatives [1]. 
Perhaps the best known one is the Riemann-Liouville def- 
inition, which is given by 
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where r(-) is the gamma function and n — 1 < a < n. 
The geometric and physical interpretation of the frac- 
tional derivatives was given in [21]. Upon consider- 
ing the initial conditions to be zero, the Laplace trans- 
form of the Riemann-Liouville fractional derivative is 
L | d | = s a L{f(t)}. So, the fractional integral op- 
erator of order "cv" can be represented by the transfer 
function F(s) — 

The standard definition of the fractional diffcrintegral 
do not allow direct implementation of the fractional oper- 
ators in time-domain simulations. An efficient method to 
circumvent this problem is to approximate the fractional 
operators by using the standard integer order operators. 
In the following simulations, we will use the approxima- 
tion method proposed in [22], which was also adopted 
in [9, 11, 14, 15, 16, 17]. In Table 1 of [9], the authors 
gave approximations for l/s q with q = 0.1 — 0.9 in steps 
0.1 with errors of approximately 2dB. We will use these 
approximations in our following simulations. 

Consider the master-slave synchronization scheme of 
two autonomous n-dimensional fractional order chaotic 
systems 
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with the master system M and the slave system S. 
Where a > is the fractional order, with which the in- 
dividual dynamical systems are chaotic, c > is the cou- 
pling strength, and T € R nxn j s a constant 0-1 matrix 
linking the coupling variables. For simplicity, we assume 
r = diag(ri, r%, ■ ■ ■ ,r„) is a diagonal matrix. If there 
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is a coupling between the zth state variable of the two 
coupled chaotic systems, then = 1; otherwise, r t = 0. 
Define the error signal as e = x — y, the aim of the syn- 
chronization scheme is to design the coupling strength 
such that ||e(t)|| — > as t — > oo. This scheme is similar 
to the master-slave synchronization of classical integer- 
order chaotic systems. 

Next, we numerically study the synchronization of frac- 
tional order chaotic systems via two examples, we first 
consider the fractional order Chua's system [9] 
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when a > 0.9, this system can produce chaotic solutions 
[9]. Particularly, when a = 0.9 and a = 12.75, the frac- 
tional order Chua's system is chaotic. The phase plot of 
x and z is shown in Fig.l. 
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FIG. 1: Phase plot of the fractional order Chua's system with 
a = 0.9. 





FIG. 2: Synchronization error of the fractional order Chua's 
systems with a = 0.9: (a) c = 4, (b) c = 7. 

We let T = diag(l,0, 0), which implies that only the 
first variable x is used to coupling the two fractional order 
chaotic systems. To obtain a critical value of c to make 
the two systems synchronized, we continuously increase 
the coupling strength c, from c = 0, in step 0.5. When 
c < 4, no synchronous phenomenon is observed. When 
c = 4, the curve of the synchronization error J(t) = 
log(||e(t)||) is shown in Fig. 2 (a), which indicate that the 
master-slave synchronization is achieved. In Fig. 2 (b), we 
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FIG. 3: Synchronization error of the integer order Chua's 
systems: (a) c = 4, (b) c = 7. 



show the curve of the synchronization error when c = 7, 
in which the synchronization effect is better than that of 
c = 4. 

For the purpose of comparison, we also plot the curves 
of synchronization error of the integer order Chua's sys- 
tems (a = 9.5) in Fig. 3. Comparing Fig. 2 with Fig. 3, 
we can know that the synchronization rate of the frac- 
tional order Chua's systems is slower than its integer or- 
der counterpart. 

We next consider the fractional order Rdssler system 
[16]: 
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when a = 0.9 and a = 0.4, the above system is chaotic. 
The phase diagram of the chaotic attractor is shown in 
Fig. 4. 

We also let T = diag( 1,0,0), and do the similar sim- 
ulations as in the above example. When the coupling 
strength c = 0.5, the two fractional Rossler systems 
achieve synchronization. The curve of the synchroniza- 
tion error of the fractional order Rossler system is shown 
in Fig. 5 (a). In Fig. 5 (b), we also plot the curve of the 
synchronization error of the integer order Rossler systems 
(a = 0.165). From Fig. 5, we know that the synchroniza- 
tion rate of the fractional order Rossler systems is also 
slightly slower than its integer counterpart. 
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FIG. 4: Phase plot of the fractional order Rossler system with 
a = 0.9. 
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FIG. 5: The curves of synchronization error: (a) the fractional 
order Rossler systems with a = 0.9 and c = 0.5; (b) the 
integer order Rossler systems with c = 0.5. 

We have also tested the synchronization scheme (2) on 
several other fractional order chaotic systems [23]. Lim- 
ited to the length of this brief report, we omit these re- 
sults here. 



In summary, in this brief report, we have studied the 
master-slave synchronization of fractional order chaotic 
systems. To our best knowledge, this is the first report 
on the synchronization of fractional order dynamical sys- 
tems. We have shown that fractional order chaotic sys- 
tems can be synchronized by utilizing the similar scheme 
as that of their integer order counterparts. 

Future works regarding this topic include the investiga- 
tion of some other types of synchronization of fractional 
order chaotic systems, such as the phase synchronization 
[24] and the projective synchronization [25], as well as the 
synchronization of fractional order hyperchaotic systems 
[16]. 
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